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Abstract 

We propose a practical way to overcome the ubiquitous problem 
of the overestimation of the zero-bias and zero-temperature conduc- 
tance, which is associated to the use of local approximations to the 
exchange-correlation functional in Density-Functional Theory when 
applied to quantum transport. This is done through partial substi- 
tution of the local exchange term in the functional by the nonlocal 
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Hartree-Fock exchange. As a non-trivial example of this effect we re- 
visit the smallest molecular bridge studied so far: a H2 molecule placed 
in between Pt nanocontacts. When applied to this system the value 
of the conductance diminishes as compared to the local-exchange-only 
value, which is in close agreement with results predicted from Time- 
Dependent Current-Density-Functional Theory. Our results issue a 
warning message on recent claims of perfect transparency of a H2 
molecule in Pt nanocontacts. 

It is widely admitted that a first-principles methodology based on Density- 
Functional Theory (DFT) is the only compromise between accuracy and 
computational resources that can help pave the way towards molecular- 
engineered nanoscale devices [1]. The DFT manageability is based, for the 
most part, in the use of the Kohn-Sham scheme (DFT-KS) [2,3]. In addition, 
the central track doing a practical DFT-KS scheme lies on the definition of 
an approach to the true exchange-correlation (xc) functional [4] . Initially, the 
manageability of DFT-KS comes upon the introduction of the Local-Density 
Approximation (LDA) to the xc functional [3]. (Henceforth, extensions be- 
yond LDA within local exchange models, like generalized gradient approxi- 
mations (GGA) [5], and Meta-GGA [6], are referred as semi-local methods 
and labeled as LDA methods for the sake of simplicity). 

This ab-initio based method combined with Landauer's formalism is nowa- 
days routinely applied to compute the conductance, G, of nanoscaled sys- 
tems, i.e metallic nanocontacts and single-molecule junctions [7-9]. However, 
the use of DFT-KS electronic structure within LDA approximations, would 
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present two problems in this regard: (i) an obvious one related to the use of 
LDA and related approximations to approach the true xc-potential, v xc [4, 9] ; 
and (ii) a much more subtle one which could be associated with the absence 
of dynamical corrections in the DFT-KS schemes [10, 11]. 

Regarding the first point, the reliability of LDA results relies naturally on 
the existence of fairly homogeneous, well-behaved electronic densities. While 
this is the case for metallic nanocontacts, molecular bridges, definitely, do 
not satisfy this premise [10]. For weakly-coupled molecules conducting in 
the Coulomb blockade regime, where charge localization is strong, failures 
of LDA are traced back to the self-interaction and derivative discontinuity 
problems, which are inherent to local or semi-local approximations [9, 12- 
14]. Interestingly, this problem is not present in more sophisticated exact 
exchange (EXX) nonlocal approaches, as recently emphasized by others [15]. 

Concerning the second point, the use of Time-Dependent Density- Functional 
Theory (TDDFT) [16] has been shown to give the exact total current if 
the exact v xc is known [10], and it reduces in the zero-frequency limit to 
the standard static treatment only when the adiabatic approximation is in- 
voked [10]. Furthermore, as initially suggested by Vignale and Kohn [11], 
the apparent impossibility to approach the exact v xc through successive local 
approximations could be circumvented by a new form of functional where the 
local current density plays a central role or, equivalently, a nonlocal density- 
based theory. Recent work [10] based on Time-Dependent Current- Density- 
Functional Theory (TDCDFT) has shown that dynamical effects manifest 
themselves in the appearance of a dynamical potential, v^ n , that always op- 
poses the electrostatic potential and lowers the conductance obtained within 
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LDA methods, even in linear response. The origin of this potential can be 
traced back to the nonlocal response to an electric field which is captured 
by TDCDFT [10, 11]. It is difficult to assess to what extent the resulting 
dynamical potential accounts for the lack of nonlocality of the starting LDA 
approximation or if it represents a true intrinsic dynamical effect which is 
present even for the exact DFT-KS potential [9] (probably both). 

The role played by the dynamical term in increasing the resistance is com- 
patible with another seemingly unrelated fact: The polarizability of molecu- 
lar systems (e.g., hydrogen [17] or polymeric [18] chains) is severely overesti- 
mated by LDA, while Hartree-Fock (HF), TDCDFT, and exact many-body 
calculations -all of them nonlocal schemes- yield similar values. We note 
here that a recent work has shown that the conductance can be expressed in 
terms of the polarizability [19]. 

In the light of the above discussion, one might ask: to what extent can 
nonlocal approaches to the v xc ( out of the realm of standard static DFT-KS 
theory) mimic the missing dynamical corrections? In this work we explore 
the influence of a nonlocal potential -such as HF exchange- on the conduc- 
tance of a molecular bridge (see Figure 1). To this end we choose to study 
electron transport in a H 2 molecule bridging Pt nanocontacts [20-22] (see in- 
set in Figure 2), using nonlocal functionals customized for this problem. We 
show that the HF-like exchange has a strong influence on the conductance 
of this bridge, rapidly decreasing as the percentage of HF-like increases in 
the hybrid functional. By adjusting this percentage, with the aid of accurate 
quantum chemistry calculations in clusters, we conclude that previously re- 
ported calculations based on LDA might have systematically overestimated 
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the value of the conductance in this system. Finally, for this benchmark sys- 
tem, we compare the numerical results derived from our method, the nonlocal 
functionals for an effective DFT-KS formalism, with those obtained within 
TDCDFT [10]. The conductance values are quite similar, which further sup- 
ports the conclusions reached here. The average value of the conductance for 
such system turns out to be rs 0.2 x 2e 2 /h = 0.2 x Go, in stark contrast to 
LDA calculations yielding w 1.0 x Go and in support of a previous work by 
some of us [21]. 

Our theoretical approach have been thus developed in accordance with 
the scheme that follows. The central track doing an effective DFT-KS scheme 
lies on the definition of the v xc [1,4,23]. Rooted on the adiabatic connection 
theorem, which formally justified the design of nonlocal functionals (usually 
referred as hybrid functionals [24]), and based on the experience gained af- 
ter original Becke's proposal [25], we adopt the following expression for an 
approach to the v xc : 

Vxc = av» F + (l-a)v l r l + v l r l , (1) 

where v^ F is the HF-like potential contribution, which accounts for nonlocal 
effects, the v 1 "™ 1 contributions are given by any of the available exchange 
(x) or correlation (c) local functionals and, finally, a is a fitting parameter. 
With this parameterized potential, we will solve the mean-field equations 
derived from DFT-KS formalism for a finite system, the cluster. This finite 
cluster is selected in such a way that the boundary conditions determined by 
the problem are fulfilled, as discussed in the next paragraphs. Finally, the 
electronic states derived from the DFT-KS are deployed in conjunction with 
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Landauer formalism for zero temperature quantum transport characteriza- 
tion [26]. Concerning the software employed for the numerical calculations, 
we use the code ALACANT (ALicante Ab initio Computation Applied to 
NanoTransport) [27] which is interfaced to GAUSSIAN03 [28]. 

We first compute the LDA conductance (a = 0) of the system shown in 
the inset of Figure 2 where the H 2 molecule lies along the Pt nanocontacts 
axis [29]. In agreement with previous works [20,22,30], we obtain a single fully 
conducting channel at the Fermi level which yields G ~ 1.0 x Go = 2e 2 /h. 
This result supports the implementation for the computation of the current 
[31]. 

To probe in a systematic way the influence of the nonlocal exchange on 
the electronic structure and concomitant change in the conductance, the 
coefficient a is now varied. Strictly speaking a should not be a constant, 
but a function of spatial coordinates, a = a(f). It is expected that HF-like 
exchange is more needed to describe electronic density regions governed by 
just one electron where the self-interaction error is more notorious [32]. This 
is particularly the case for the H 2 molecule. The value of a, on the other hand, 
should decrease as we move into the bulk electrodes where LDA or GGA is 
expected to perform better. As a compromise between both situations we 
take a to be constant in the nanocontacts region close to the H 2 molecule 
where the coordination is low, i.e., the bridge Pt-H 2 -Pt (see inset in Figure 
3). The rest of the system is described by the self-energy corresponding 
to a parameterized Bethe lattice. Figure 1 shows the effect of the HF-like 
exchange on the conductance. For values close to 50% the conductance has 
almost dropped to zero. 
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The drop in conductance as the amount of HF-like exchange increases 
at the bridge can also be understood from a purely microscopic point of 
view. This construction is based on the approach to the exact non-correlated 
density-based scheme for determine the electronic structure [23], and its com- 
bination with the accepted Landauer formalism for determine the electronic 
conductance [7]. Therefore, under these two conditions, the electronic struc- 
ture limited to the HOMO-LUMO gap, A H omo-lumo, will play a key role in 
the final outcome of the electronic conductance [33] . Figure 3 shows the effect 
of the HF-like contribution on the HOMO-LUMO gap (A H omo-lumo) for 
the isolated Pt-H 2 -Pt bridge. This effect is well known in molecules and bulk 
semiconductors where LDA approximations usually underestimate the real 
value of the charge gap [34-36]. The one electron Green's function (g(e)) 
of the isolated Pt-H 2 -Pt bridge has poles at the position of the molecular 
orbitals energies. Once the isolated system is coupled to the rest of the elec- 
trodes these poles shift and broaden in a way dictated by the self-energies 
El and Er representing the electrodes: 



The Fermi level is likely to lie in the gap between the HOMO and LUMO 
states of the Pt-H 2 -Pt cluster since no significant charge transfer is expected 
between Pt atoms. As can be seen in Figure 3, the LDA gap is almost zero 
for this bridge. From the opening of the gap as a increases one could easily 
have anticipated a strong decrease in conductance since this depends on the 
Green's functions through the well-known expression [26]: 
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W(e) = [(c±t < 5)/- J ff-E L -E R ]- 1 . 



(2) 



G 



2e2 



Tr[ r L (e) g+(e) T R (e) <T(e)]. 



(3) 
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Finally, to estimate a we follow a well-established methodology in semi- 
conductors and insulators. The value is chosen so that Ahomo-lumo matches 
the charge gap, GAP. Accordingly, it is defined as GAP = 2E {N)-E (N + 
1) — E (N — 1), where N represents the number of electrons, and E is the 
total energy for the selected ground states in a Kohn-Sham system of indepen- 
dent electrons. Matching between Ahomo-lumo and GAP will be enough 
to guarantee accuracy for mean-field Green's functions [37,38], and therefore 
for the conductance value, expressions 2 and 3. The choice of the finite region 
where to perform the search of a is constrained by three conditions: (i) it 
must contain the zone which is restricting the current flow, in our case the H 2 
molecule; (ii) the HOMO and LUMO must not be localized on the electrodes; 
and (iii) screening effects from the contacts should be -at least- partially in- 
cluded. These three conditions reduce the range of possibilities to the cluster 
shown in Figure 3, and limit the selection of system (previously referred as 
the extended molecule [9,39]). One way to estimate to what extent screening 
is included in such small cluster is to look for a quantity such as the internal 
mode frequency of the H 2 molecule, wh 2 , and check how much it changes 
with increasing the size of the cluster. We find that larger clusters do not 
appreciably change that quantity, which is given by wh 2 ~ w h 2 /^ where w H2 
is the internal mode of an isolated H2. This change in the frequency with 
respect to the free molecule suggests an intermediate regime of coupling. 

To find out ground state energies and therefore the right GAP, we begin 
by considering the mean-field HF method. To approach the exact many-body 
solution to the GAP, we invoke hierarchical correlated ab initio methods 
such as Moller-Plesset perturbation (MP) or Coupled-Cluster (CC) theories. 
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The corresponding results are presented in Figure 4. The GAP for this 
system, ~ 5 eV, critically differs from the Ahomo-lumo, described from 
the electronic structure in pure DFT-LDA calculations ~ eV, see Figure 3. 
Therefore, and considering the scheme we propose to correct this pitfall, the 
Ahomo-lumo for the system is extracted from Figure 4. Comparing those 
results with the whole information contained in curve of Figure 3, we estimate 
that the percentage of HF exchange in Eq. 1 should be ~ 60 percent, which 
gives a very low value for the conductance (see Figure 1). It is the main 
result we should draw from this paper. 

Note that there is an increasing concern about the influence of exact 
exchange for conduction phenomena at the molecular level [40,41]. Although 
a rule of thumb for the needed admixture of exact exchange is provided here, 
we indeed believe that it will be difficult to unambiguously give the optimum 
weight of exact exchange for an uniformly good description of all systems. 
Instead of this, we would like to remark how the results clearly show that 
inclusion of exact exchange is necessary for accurate descriptions. We should 
point that the same amount of HF exchange is obtained (with error band of 
10 percent), if the contour condition to find out its value is limited to the 
matching between the HOMO eigenvalue and the negative of the ionization 
potential, enoMO = —IP [42]. Concerning future work and despite being 
technically challenging, it would be definitively interesting to investigate the 
performance on these benchmark systems of a new generation of functionals 
describing v xc . For instance, we should mention orbital-based functionals [43] , 
self-interaction corrected functionals [44,45], the hybrid functionals based on 
a screened Coulomb potential [46,47], and the last generation of v xc density- 
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functional belonging to the M06 suite of methods [48,49]. Furthermore, 
more studies for closely related molecules and nanocontacts will be highly 
welcomed too. 

If we were able to follow unambiguously the same procedure for larger 
systems the value for a would probably be reduced, although we do not 
expect it to go to zero. Unable to get rid of this uncertainty we compare our 
results to those obtained with TDCDFT following the work by Sai et al. [10]. 
A dynamical counteracting potential appears whenever there is a change in 
the macroscopic density which modifies the local-density conductance: 

n — ^local /.\ 

1 + G local R d y- y } 

From the average density profile along the z-direction (see Figure 5) we 
compute the corrections to the local conductance. The electronic density 
is fixed in the electrodes to the bulk value, considering for the number of 
free-electrons as the mean experimental atomic-point-contact transmission. 
In case of Pt-nanocontacts, this number implies a range between one- and 
two-free-electrons per Pt atom [20]. The rest of the parameters involved in 
the calculation of R dyn were estimated in agreement with the method pro- 
posed by Sai et al. [10]. Our calculations brings the conductance down to 
~ 0.2 x Go- (Details are included in Ref. [50]). Nevertheless, this procedure 
is not free from ambiguities either, but it points towards similar conclusions 
as the ones drawn above. 

In summary, we have presented a route to approach the DFT-KS elec- 
tronic transport calculations based on the controlled addition of EXX con- 
tribution to v xc . The selection of the fraction of HF-like needed is based on 
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first-principles procedures. Our results agree with the ones obtained with the 
method developed by Sai et al [10] within TDCDFT, and supports the func- 
tionality of DFT methods based on nonlocal models for v xc . Its agreement 
opens up possibilities to formally justify a connection between nonlocal DFT 
and TDCDFT in the zero-frequency limit. Furthermore, more work is needed 
on this respect since the exact amount of nonlocality might be difficult to 
estimate given: (i) the difficulties for the definition of the region where it 
should operate; and (ii) the correct implementation of this behavior. How- 
ever, we can at least define an upper limit for this nonlocal HF correction, and 
give arguments in favour of its non-zero value, as DFT-LDA methods consid- 
ered. Within this context, the numerical calculation of the conductance for 
a H 2 molecule in Pt nanocontacts have been corrected, lowering in one order 
of magnitude the standard DFT-LDA calculation. Note, however, that the 
statements and the conclusions that will be made are not limited to calcu- 
lations of conductance in this benchmark system. It is widely accepted that 
nonlocal contributions to v xc do play a significant role to get beyond DFT- 
LDA results either in molecular systems or bulk materials. For instance, 
we could mention the polarizability of polymeric chains [18], the magnetic 
coupling and the band structure in systems such as antiferromagnetic in- 
sulators [51,52] and magnetic oxides [53]. Finally, our theory is that the 
over estimation of the conductance is not intrinsic to the DFT-KS approach. 
This wrong behavior could be greatly improved by a controlled design of an 
exchange-correlation potential which incorporates nonlocal contributions. 
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Figure 1. Reduction of the conductance due to the progressive sub- 
stitution of the local xc-potential by the nonlocal HF-like contribution. 
The inset shows the conductance as function of energy calculated in a 
nonlocal DFT approximation for a = 0.6. 

Figure 2. Transmission as a function of energy calculated in a local 
DFT approximation for the Pt-H 2 -Pt bridge shown in the inset. 

Figure 3. Enhancement of the HOMO-LUMO gap {^homo-lumo) 
due to the substitution of local xc potential by the nonlocal HF con- 
tribution. The inset shows the Pt-H 2 -Pt system bridging the Platinum 
electrodes. 

Figure 4. Ground state energy calculations of the charge gap, GAP, 
for the system shown in the inset. GAP = 2E (N)-E (N+1)-E (N- 
1), where N represents the number of electrons, and E is the total en- 
ergy for the selected ground states. Total energy calculations were 
done using wave- function based methods: Hartree-Fock (HF), Moller- 
Plesset perturbation theory at the second, third and fourth order (MP2, 
MP3, MP4, respectively) and Coupled-Cluster theory with single, dou- 
bles and single, doubles, and perturbatively estimated triples [CCSD, 
CCSD(T), respectively] methods. 

Figure 5. Electron density profile after planar average (dashed line) 
and three-dimensional average of the microscopic density. 
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Figure 1 
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Figure 2 
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